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Universality and Lepton Flavor Violation in T Decays at CLEO 

J.E. Duboscq a , for the CLEO Collaboration * 

a Wilson Laboratory, Cornell University, Ithaca, NY, 14853, USA 

I present two analyses done with the CLEOIII detector at CESR. The first is a search for the LFV decay 
T — > fir with preliminary results, and the second constrains Lepton Universality in T — > rr relative to T — > n/i. 
This second analysis, whose results are final, has dramatically improved uncertainties relative to other such 
measurements, and is also the first observation of the decay T(3S) — > rr. A limit is also set on the involvement 
of a CP odd Higgs in the decay of the T(IS') to r pairs. 



1. Introduction 

In this report, I present results on two analyses 
done with data taken using the CLEOIII detector 
on the CESR storage ring by the CLEO Collabo- 
ration. The first analysis is a search for the Lep- 
ton Flavor Violating (LFV) decay T(nS) — > fir, 
n = 1, 2, 3 - all results from this analysis are pre- 
liminary. The second analysis tests Lepton Uni- 
versality in the decays T — > (j,[i,tt, n = 1,2,3 - 
all results from this analysis are final. 

The CLEOIII detector is a modern almost her- 
metic detector, sitting on the CESR beamline, 
where positrons and electrons collide with an en- 
ergy « 5 GeV. The detector included a silicon 
strip detector at its center, surrounded by a wire 
chamber, which provided tracking and dE/dx in- 
formation, a Cesium Iodide calorimeter, and a 
Ring Imaging Cerenkov Detector all in a 1.5 Tesla 
magnetic field. Outside of this were a series of 
muon chambers. Further details of the detector 
can be found in [Q. 

The data used in this analysis comprise the 
T(nS), n = 1, 2, 3 datasets, along with a portion 
of the T(4S) data for validation purposes. These 
datasets include data taken at the peaks of the 
resonances, as well as smaller datasets taken some 
30 — 40MeV below the resonances. The total 
sample includes approximately 20 x 10 6 T(15), 
10 x 10 6 T(25) and 5 x 10 6 T(35) decays. 



2. Search for LFV in T — * fir 

Lepton Flavor Violation (LFV) might provide 
a key to understanding Lepton Number Viola- 
tion, as well as Baryon Number Violation, and 
hence might be helpful in understanding the 
matter/anti- matter asymmetry in the Universe. 
LFV decays in the T sector can be related to LFV 
decays of the r lepton by a simple reordering of in- 
put and output lines in the Feynman diagram for 



the process e H 



r* 



TfJ,. 



In this way, 
3fx) < 10~ 6 



|2j have show that a result of B(t 
implies that B(T -> (xt) < 10~ 2 . 

LFV decays of the T to lepton pairs might also 
be expected in the presence of SUSY loops [3], 
or other exotic physics such as leptoquarks and 
theories with more than one Z boson. 

A general ansatz for talking about LFV decays 
of the T involves a generic four fermion interac- 
tion vertex between two b quarks, a muon and a 
r. With the coupling constant at the interaction 
denoted a at and the relevant new physics mass 
scale denoted A, one finds 0]: 



B(T -> tt) 



oc (a N /a) 2 (M T /A) 4 



*I would like to thank A.Lusiani and the staff of Tau06 
for a wonderfully stimulating conference. 



The CLEO analysis presented here is a search 
for T(nS) — > fir, r — > evv, n=l,2,3. The final 
observed state is a two track event with a muon 
and a electron and missing energy. The muon 
would have an energy a little below the beam en- 
ergy, while the electron energy spectrum would 
be approximately that expected for electrons in 
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normal r pair events. 

The fit to the data is performed using an ex- 
tended maximum likelihood function, composed 
of a product PDF, summed over expect signal 
LFV shapes, direct ee — > rr shapes, ee — > 
and ee — > \x\x with — > ew decays. Data from 
the T(4S I ) resonance, and off resonance data are 
used as calibration and control samples. 

Fig n shows the beam energy scaled electron 
momentum spectrum versus the scaled muon mo- 
mentum spectrum, in T(45) data where no signal 
is expected, since even T(4S I ) — * rr is too small 
to be observed. 




Figure 1. Beam energy scaled electron ver- 
sus the scaled muon momentum in T(AS) data. 
The left hand side is dominated by r pair de- 
cays. The signal region is a vertical strip near 
Pfj,/ Eseam ~ 0.97. The concentration of events 
near p^/EBeam = 1 and p e /E Beam = 0.5 is from 
/Z/X7 backgrounds. Also at p^/EBeam — 1 are 
events from muon decays to electrons in pro- 
duction. 



Fig [21 shows the data distributions for the fit 



quantities at the T(15). Note the clear absence 
of a signal, especially in the muon momentum 
spectrum. 



* 




Figure 2. Likelihood distributions for the scaled 
muon momentum, E/p, dE/dx and scaled mo- 
mentum for the electron candidate. Note the 
curve indicating what 100 signal events would 
look like. 



Preliminary analysis results are summarized in 
Tabled The largest systematic uncertainties are 
from PDF shapes and their correlations. This 
study represents the first limits on B(T — > r/i). 
These limits set a lower limit of « 1 TeV on the 
scale of new physics, assuming a strong coupling. 

3. Lepton Universality in T — > /u/i, rr 

In the Standard Model, the couplings between 
leptons and gauge bosons are independent of the 
lepton flavor, so the branching fractions for the 
decay T(nS) — * should be independent of 

the flavor of the lepton I, except for negligible fi- 
nal state lepton mass effects. Any deviation from 
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Resonance 


T(IS) 


T(25) 


T(35) 


Efficiency 


8.9% 


8.9% 


8.9% 


Events 


< 10.0 


< 10.7 


< 8.5 


B{T -> T fi) (10- b ) 


< 6.2 


< 25 


< 22 


B{tu)/B{m) 


< 0.023% 


< 0.17% 


< 0.13% 



Table 1 

Preliminary results for the LFV search. Displayed are the signal MC efficiencies, raw event counts, 
resulting branching fractions, and the ratio B(T — > t/i)/B(T — > /i/i). All limits are 90% confidence level 
upper limits. 



unity for the ratio of branching fractions 7Z^ T = 
B(T(nS) — ► tt) / B(T(nS) — ► /i/i) would indicate 
the presence of new physics. The ratio 1Z^ T is sen- 
sitive to the mechanism proposed in jo] , in which 
a low mass CP-odd Higgs boson, A , mediates 
the decay chain T(15) — * i]bj, r)b — > A — * tt. 

CLEO has recently measured the partial width, 
r ee , from e + e~ — > Y(nS') 6 , n = 1,2,3, as well 
as the branching fraction for T(nS) — > /i + M~ 13- 
This analysis complements these measurements 
by measuring 7^ T directly, and scales this re- 
sult to obtain B(T(nS) — > tt). An upper limit 
on the product branching fraction B(Y(IS) — > 
r]b^)B(r]t — > A — > tt) is extracted. 

The analysis technique, similar to that in [7j, 
isolates the T — * /i/i, tt signals by subtracting 
a luminosity and beam energy weighted number 
of events observed in off-resonance data from the 
number observed in on-resonance data, and, after 
further background correction, attributes the re- 
maining signal to T decays to leptons. Selection 
criteria are developed to isolate /i/i and tt final 
states using a subset of the data acquired near 
the T(4S). Another subset of T(45) data is used 
to verify that subtracting the scaled off-resonance 
data from the on-resonance data produces no sig- 
nal for Y(4S) — * tt, /i/i, indicating that non-T 
backgrounds are suppressed by the subtraction. 
As a further crosscheck, the off-resonance produc- 
tion cross-sections for tt and \x[L are verified to 
agree with theoretical expectations. 

The final states chosen for both the T — ► /i/i 
and T — > tt decays are required to have ex- 
actly two good quality charged tracks of opposite 
charge. 



Selection criteria for the final state closely 
parallel those of 7 , requiring tracks with mo- 
menta scaled to Bbcam between 0.7 and 1.15, of 
which at least one is positively identified as a 
muon. The energy deposited by a particle in a 
calorimeter shower, Eqc, is required to satisfy 
100 MeV < E C c < 600 MeV. No more than one 
shower unassociated with a track and with energy 
above 1% of -Ebcam is allowed. 

At least two neutrinos from final states of 
e + e~ — > tt escape detection. Thanks to CLEO- 
III's hermeticity, the following criteria select such 
events despite the energy carried away by the un- 
reconstructed neutrinos. The total charged track 
momentum transverse to the beam direction must 
be greater than 10% of i?bcam, and the total 
charged track momentum must point into the 
barrel region of the detector where tracking and 
calorimetry are optimal. Events with collinear 
tracks are eliminated. Tracks are required to have 
momenta greater than 10% of -Ebeam to ensure 
that they are well-reconstructed, and, to mini- 
mize pollution from two-particle final states, they 
are required to have momenta less than 90% of 
-Ebeam- The total observed energy due to charged 
and neutral particles in the calorimeter is simi- 
larly required to be between 20% and 90% of the 
total center-of-mass energy. To reduce overlap 
confusion between neutral and charged particles, 
a shower's energy scaled to its associated track's 
momentum must be less than 1.1. 

Final states are further exclusively divided ac- 
cording to the results of particle identification 
into (e,e), 0,e), (e,^), 0,/i), (e,A), (fi,X), 
(X,X) sub-samples, with particles listed in de- 
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scending momentum order. The first (second) 
particle listed is referred to as the tag (signal). 
The sub-samples are used to crosscheck consis- 
tency of results across decay modes. Lepton 
identification requires a track momentum greater 
than 500 MeV to ensure that the track intersects 
the calorimeter. Electrons are identified by re- 
quiring that 0.85 < Ecc/P < 1-10, where P is 
the track momentum, and that the specific ion- 
ization along the track's path in the drift chamber 
be consistent with the expectation for an elec- 
tron. A muon candidate in r decays is a charged 
track which is not identified as an electron, having 
momentum above 2 GeV (1.5 GeV) for a tag (sig- 
nal) track and confined to the central barrel where 
beam related background is a minimum. Further- 
more, the energy deposited in the calorimeter for 
this track must be between 100 and 600 MeV, and 
the particle must penetrate at least three inter- 
action lengths into the muon detector. Particles 
identified as neither e nor \i are designated X, and 
are a mixture of hadrons and unidentified leptons. 

The decay products of r pairs from T decays 
tend to be separated into distinct hemispheres. 
Since the photon spectrum expected in r decays 
depends on the identity of the charged particle, 
calorimeter showers are assigned to either the tag 
or signal hemisphere according to their proximity 
to the tag side track direction. For each r decay 
mode pair, the number of unmatched showers in 
the calorimeter as well as the total energy of these 
showers on each side of the decay are used as se- 
lection criteria. 

For the (e, e) and (//, /i) modes, the sum of the 
magnitudes of the tag and signal track momenta 
must be less than 1.5i?beam, reducing the con- 
tamination from radiative dilepton events. To 
reduce backgrounds from e + e~ — > and 
e + e~ — > e + e~l + l~ in the (e,e), (/z, fi), [X, X) 
categories, the minimum polar angle of any un- 
seen particles, deduced from energy-momentum 
conservation, is required to point into the barrel 
region, where calorimetry cuts will ensure rejec- 
tion. 

Potential backgrounds due to cosmic rays are 
accounted for as in [Jj. In all cases these were 
negligible. 

Figure shows the superimposed on-resonance 
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Figure 3. Total energy distributions, scaled to 
center of mass energy, E TT /^/s for the tt final 
states at the Y(nS), n= 1,2,3,4 on (points) and 
scaled off (histograms) resonance. The excess of 
on-resonance relative to scaled off-resonance data 
is attributed to T decays for the lower resonances, 
while the agreement at the Y(45) tests the valid- 
ity of the subtraction. 



and scaled off-resonance total energy distribu- 
tions for the tt sample for all resonances. The 
scale factor is S = {C 0n / £off)(soff / son)<5intorf , 
where C and s are the data luminosity and 
squared center of mass collision energies on and 
off the resonances, and 5i n terf is an interference 
correction. The luminosity is derived from the 
process e + e~ — > 77 |5], which does not suffer 
backgrounds from direct T decays. The inter- 
ference correction <5i n tcrf accounts for the small 
interference between the process e + e~ — ► 11 and 
e + e~ — * T — > II and is estimated to be 0.984 
(0.961, 0.982) at the T(15) (25, 35) and negli- 
gible for the T(45). Note that the interference 
largely cancels in the ratios considered in this 
work. The agreement of the distributions for the 
T(45) validates the subtraction technique, and 
also highlights the absence of any process whose 
cross-section does not vary as 1/s. This agree- 
ment extends to the individual sub-samples. 

The ratio of measured relative lepton pair pro- 
duction cross-sections, TZ^ — oVt/cV/j, with re- 
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spect to that theoretically expected at the off- 
resonance points, 72° r ffthc ° ry = CT^'y/o-theory 

is K°?/1l°? thcoiy = 0.96 ± 0.03 (0.97 ± 0.03, 
0.97±0.03 , 1.00±0.03) below the T(1S) (25, 3S, 
<tS) for the sum of all t decay mode pairs, with 
statistical and systematic uncertainties added in 
quadrature. The expectation T^Off theory = .83± 
0.02 (syst) 2 , which is found to be numerically in- 
dependent of the particular resonance considered, 
and the reconstruction efficiencies are derived 
from the FPair [Hj and Koralb/Tauola |llll2H3j 
Monte Carlos. Backgrounds were corrected by us- 
ing e + e~ — > qq(q = u, d, c, s) Monte Carlo simula- 
tions HEI [13 IH- The scatter in the central 
values of ftOff/^Off theory indicates that system- 
atic uncertainties are small. 

The reconstruction efficiency for observing 
T -> fifi is derived from the CLEO GEANT- 
based simulation JSj JS] ^21 |14| . as shown in 
Table [21 This efficiency is found to be constant 
across the resonances. 

The reconstruction efficiency for observing 
T — > tt is derived using the Koralb/Tauola event 
generator integrated into the Monte Carlo simu- 
lation. Although this generator models the pro- 
cess e + e~ — > 7* — > t + t~ , the quantum numbers 
of the T and 7 are the same so it can be used 
as long as initial state radiation (ISR) effects are 
not included. This efficiency is also found to be 
consistent across all resonances within any given 
tt decay channel. 

Results of the subtraction are summarized in 
Tabled showing the first observation of T(3S) — > 

TT. 

Backgrounds resulting from cascade decays 
within the bb system to 11 are estimated using 
the Monte Carlo simulation, with branching frac- 
tions scaled to the values measured in this study. 
Cascade backgrounds with non and non tt fi- 
nal states are estimated directly from the Monte 
Carlo simulation. 

Figure |H displays the off- resonance subtracted 
data, superimposed on Monte Carlo expectations. 
The distributions shown are a sampling of r de- 



2 This expectation is lower than 1 because of the larger 
phase space available for initial state radiation production 
of fifj, relative to tt final states. 



cay modes for the momenta of the signal and tag 
tracks, as well as the total reconstructed energy. 
In all cases the Monte Carlo expectations are con- 
sistent with the data assuming lepton universality 
and branching fractions as measured in [7]. The 
agreement across the various kinematic quantities 
indicates that backgrounds are well controlled. 

Figure0shows the agreement across all tt sub- 
samples of the ratio of off-resonance cross sections 
for tt and /i/i production, relative to expecta- 
tion, as well as the ratio of branching fractions for 
each of these decay modes at the different T res- 
onances. The agreement across tt sub-samples 
both on and off the resonances is again an indi- 
cation that backgrounds are small and well esti- 
mated. 

The ratio of branching fractions and final 
branching fractions are listed in Table [31 These 
results show that lepton universality is respected 
in T decay within the rs 10% measurement un- 
certainties. 

Systematic uncertainties, summarized in Ta- 
ble 01 are estimated for the ratio of branching 
fractions, and for the absolute branching fraction. 
The ratio of branching fractions is insensitive to 
some common systematic uncertainties. For in- 
stance, the uncertainty on lZj r due to a conser- 
vative 1% variation in the scale factor S, as de- 
termined using the process e + e~ — > near the 
T(45) resonance, is found to be 0.4% or less. 

Most systematic uncertainties due specifically 
to 11 selection are derived by a variation of the 
selection criteria over reasonable ranges in the 
T(15) sample, which has the lowest energy re- 
leased in its decay. The most significant of these 
are due to momentum selection (1.3%), calorime- 
ter energy selection (1.1%) and angular selection 
(1.1%). The systematic uncertainty due to mod- 
eling of the trigger, also included in the 11 se- 
lection, is estimated by using a loose pre-scaled 
tracking trigger and comparing it to the more so- 
phisticated triggers used in this analysis. This 
variation leads to a systematic uncertainty esti- 
mate of 1.6%. 

Backgrounds are assumed to be due solely to 
T decays. As in [Jj, these were chiefly due to 
cascade decays to lower resonances, and are esti- 
mated to be 2.5% (15%, 11%) of the tt sample at 
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T(15) 


T(2S) 


T(35) 


7V(mm) (10 3 ) 


345 ±7 


121 ±7 


82 ±7 


W (%) 


65.4 ± 1.2 


65.0 ±1.1 


65.1 ±1.2 


JV(T - MM) (10 3 ) 


527 ± 15 


185 ±11 


126 ± 11 


N(tt) (10 3 ) 


60.1 ±1.5 


21.8 ± 1.5 


14.8 ± 1.5 


e(rr) (%) 


11.2 ±0.1 


11.3 ±0.1 


11.1 ± 0.1 


N(T -> tt) (10 3 ) 


537 ± 14 


193 ±12 


132 ± 13 



Table 2 

Summary of reconstructed events for T — > MM (top) and T — > tt (bottom). Shown are the number of 
events (Nit) after subtraction of backgrounds estimated from scaled off-resonance data and T feed-through 
estimated from the Monte Carlo simulation, the signal efficiency (e(ll)), and the total efficiency corrected 
number of signal events N(T — > 11) = Nu/e(ll). The tt events are summed over all decay modes of the r. 
Uncertainties included in this table include data and Monte Carlo statistical uncertainties, uncertainties 
on backgrounds, and detector modeling (included only for e(mm) to avoid double counting in the final 
ratio) . 





;v Tr 


B(T -> TT ) (%) 


T(1S) 


1.02 ±0.02 ±0.05 


2.54 ±0.04 ±0.12 


T(25) 


1.04 ±0.04 ±0.05 


2.11 ±0.07 ±0.13 


T(3S) 


1.07 ±0.08 ±0.05 


2.55 ±0.19 ±0.15 



Table 3 

Final results on the ratio of branching fractions to tt and MM final states, and the absolute branching 
fraction for T — > tt. Included are both statistical and systematic uncertainties, as detailed in the text. 



the T(1S)(2S, 3S 1 ), with an estimated uncertainty 
contribution to TlJ T of 0.1% (2.4%, 1.3%). 

The uncertainty due to detector modeling in jj] 
was estimated to be 1.7%: this value is used here 
conservatively for the systematic uncertainty on 
the ratio. 

The modeling of the physics in T(IS') — > fifj,, 
obtained by varying the decay model for T — > MM 
between the Monte Carlo simulation and Koralb 
with ISR simulation turned off, contributes an un- 
certainty 2%. This uncertainty is consistent with 
the variation in the product of the off-resonance 
cross section and reconstruction efficiency using 
the FPair, Koralb, and Babayaga Monte Carlo 
simulations, and is thus likely conservative, as di- 
rect MM production from the T at the peak in- 
volves much lower energy final state photons than 
off resonance production. An uncorrelated un- 
certainty of 2% for modeling of T — > tt is as- 



sumed, consistent with the uncertainty on the 
off-resonance production cross section derived in 
previous analyses, and is again conservative as on- 
resonance production of tt final states involves 
fewer photons than direct continuum production. 
To test the sensitivity to ISR simulation, the re- 
construction efficiency for events with no ISR sim- 
ulation is compared to that for events generated 
with ISR simulation turned on and re-weighted 
according to the relative value of the T line shape 
at the t pair mass. These two efficiencies agree 
to within 0.8% of their central value. 

The mechanism described in [5] could induce 
a value of TZ^ T not equal to one. By assuming 
that the mass of the 776(15*) is 100MeV/c 2 be- 
low the T(IS') mass, consistent with the largest 
value in the value quoted for TZ^ T (\S) can be 
translated into an upper limit on the combined 
branching fraction of B(T(1S) — > f?&7)-B(?7b — ► 
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Figure 5. Breakdown by mode of off- and on- 
resonance data at the different resonances. On 
the left, the ratio of the production cross sec- 
tion for e + e~ — > = t,/j,), relative to its ex- 
pectation, is plotted for data taken below the 
T for each r decay mode pair. On the right, 
the ratio of branching fractions for the process 
T — * = r, fi) is displayed. The line centered 
at 1 in each case represents the Standard Model 
expectation. Errors shown are statistical. 



Figure 4. Distributions for the rr final states 
at the T(nS),n = 1,2,3 after subtraction of S- 
scaled off-resonance data. Distribution a) shows 
Psig/Ebeam m T(15) decays for the sum of r de- 
cay modes including exactly two identified lep- 
tons. Distribution b) shows Ptag/^bcam m T(2S) 
decays for the sum of r decay modes including ex- 
actly one identified lepton. Distribution c) shows 
E TT /y/s for T(35) for the sum of all r decay 
modes. In all cases, the solid line shows the ex- 
pected total signal and background distributions, 
assuming lepton universality. In distribution c), 
the signal and total background distributions are 
explicitly displayed. Data uncertainties shown 
are purely statistical. 



A -► rr) < 0.27% at 95% confidence level. 
Since the transition photon is not explicitly re- 
constructed, this limit is valid for all r)b that ap- 
proximately satisfy M(T(1S)) -M(n b )+T(r/ b )) < 
O(100 MeV/c 2 ). 

In summary, using the full sample of on- 
resonance T(nS), n = 1,2,3, collected at the 
CLEO-III detector, we have made the first obser- 
vation of the decay T(35) — » rr. We have also 
reported the ratio of branching fractions of T de- 
cays to rr and final states, and find these to be 
consistent with expectations from the Standard 
Model. These ratios have been combined with re- 
sults from [xj to provide absolute branching frac- 
tions for the process T — * rr, shown in Tabic 
resulting in the most precise single measurement 
of B(T(1S) — ► rr) [T^], a much improved value 
of B(T(2S) — > rr) and a first measurement of 
B(T(3S) — > rr). The ratio of branching fractions 
for rr and ^/x final states has also been used to 
set a limit on a possible Higgs mediated decay 
window. 
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Source 


Csj/st (%) 


S 


0.2 / 0.4 / 0.3 


Background 


0.1 / 2.4 / 1.3 


t, fi Selection 


2.9 


T -> n/j, Model 


2.0 


T -> tt Model 


2.0 


Detector Model 


1.7 


MC Statistics 


1.9 / 1.0 / 1.0 




4.8 / 4.4 / 4.6 


a(B TT )/B TT 


4.0 / 3.8 / 3.9 



Table 4 

Summary of systematic uncertainties. The entry 
(t(TZ^ t )/TZ^ t indicates the relative uncertainty on 
TZ^t' while the a(B TT )/ B TT entry indicates uncer- 
tainties specific to r decay modes used in addition 
to those in [7j to obtain B(T — * tt). The uncer- 
tainty on S and the background are included in 
the statistical uncertainty only. Lines with three 
entries indicate the contribution from the Y(IS'), 
T(2S) and T(3S), respectively. 
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